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We present a class of exarf cosmological solutions of the low energy string effective action in 
the presence of a homogeneous magnetic fields. We discuss the physical properties of the obtained 
(fully anisotropic) cosmologies paying particular attention to their vacuum limit and to the possible 
isotropization mechanisms. We argue that quadratic curvature corrections are able to isotropize fully 
anisotropic solutions whose scale factors describe accelerated expansion. Moreover, the degree of 
isotropization grows with the duration of the string phase. We follow the fate of the shear parameter 
in a decelerated phase where, dilaton, magnetic fields and radiation fluid are simultaneously present. 
In the absence of any magnetic field a long string phase immediately followed by radiation is able to 
erase large anisotropies. Conversely, if a short string phase is followed by a long dilaton dominated 
phase the anisotropies can be present, in principle, also at later times. The presence of magnetic 
seeds after the end of the string phase can induce further anisotropies which can be studied within 
the formalism reported in this paper. 



I. FORMULATION OF THE PROBLEM 

Suppose that at some time t\ the Universe becomes transparent to radiation and suppose that, at the same time, 
the four-dimensional background geometry (which we assume, for simplicity, spatially flat) has a very tiny amount of 
anisotropy in the scale factors associated with different spatial directions, namely 



ds 2 = dt 2 - a 2 (t)dx 2 - b 2 (t)[dy 2 + dz 2 



(1.1) 



where b(t), as it will be clear in a moment, has to be only slightly different from a(t). The electromagnetic radiation 
propagating along the x and y axes will have a different temperature, namely jl],|^ 



rri /i\ rp °1 rp — f H(t)dt rp ( ,\ rp ^1 rp - f F(t)dt 

T x (t) = Tx — = Tie J y ' , Ty{t) = Ti— = T x e J , 
a o 



(1.2) 



where H{t) and F(t) are the two Hubble factors associated, respectively with a(t) and b(t). The temperature 
anisotropy associated with the electromagnetic radiation propagating along two directions with different expansion 
rates can be roughly estimated, in the limit H(t) — F(t) <C 1, as 



^ ~ J [H(t) - F(t)]dt = lj e(t)dlogt 



(1.3) 



where, in the second equality, we assumed that the deviations from the radiation dominated expansion can be written 
as F(t) ~ (1 — e(t))/2t with |e(f)| <C 1. As we will see in the following e can be connected with the shear parameter. 

The dynamical origin of the primordial anisotropy in the expansion Q can be connected with the existence of 
a primordial magnetic field (not dynamically generated but postulated as an initial condition) or with some other 
sources of anisotropic pressure and, therefore, the possible bounds on the temperature anisotropy can be translated 
into bounds on the time evolution of the anisotropic scale factors. 

In the context of string cosmology Q there are no reasons why the Universe should not be mildly anisotropic as 
it was recently pointed out g|. Indeed, the Universe can be anisotropic as a result of the Kasner-like nature of the 
pre-big-bang solutions. Moreover, we could argue, that the Universe has to be only mildly anisotropic. By mildly 
anisotropic we mean that the scale factors should be both expanding (in the String frame). The low energy beta 
functions can be solved for a metric (1.1) with the result || 



o(t) 



t ' 


a. 


t ' 


7~i. 


, b(t) = 





<j>{t) = (a + 2(3 - 1) log 



(1.4) 
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where o? + 2/3 2 = 1. If we want a > and b > (expanding solutions) and a > 0, b > (inflationary solutions) the 
exponents a and /? have to lie in the third quadrant along the arc of the ellipse reported in Fig. [I]. 
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FIG. 1. The axionic logarithmic energy spectra (in frequency) are reported for different dilaton-driven models with Kas- 
ner-like exponents (a, 13) in the string frame. The models belonging to the third quadrant and localized on the arc of the 
"vacuum" ellipse a +2/3 = 1 correspond to anisotropic models with expanding and accelerated scale factors. If a and f3 lie 
on the ellipse but in the first quadrant we have solutions of the low energy beta functions which are both contracting (in the 
String frame) . If a and on the ellipse but either in the second or in the fourth quadrant we have solutions of the low energy 
beta functions where one of the two scale factors expands and the other contracts. Notice that the dashed (thick) line does 
correspond to the case of fully isotropic solutions (i.e. a = f3 = —l/s/S) whose intersection with the vacuum ellipse lies in 
the region of red spectra. Above the dot-dashed (tick) line the dilaton decreases for t < whereas below the dot-dashed line 
the dilaton increases. We would be tempted to speculate that to have an increasing dilaton is a sufficient condition in order 
to have a pre-big-bang dynamics. This is not the case. Indeed, an increasing dilaton is also compatible with a and /3 in the 
second quadrant where the scale factors are not both expanding. Therefore, an increasing dilaton is not a sufficient condition 
for anisotropic pre-big-bang dynamics but only a necessary condition. 



The requirement of having, in the same class of fully anisotropic models, axion spectra decreasing at large distances 
(i.e. blue frequency spectra) will select a slice exactly in the same region of parameter space ||. Therefore, general 
considerations seem to point towards mildly anisotropic pre-big-bang models. Since for phenomenological considera- 
tions H we would like to have blue (or flat) logarithmic energy spectra it seems reasonable to analyze fully anisotropic 
string cosmological models. Recently it was pointed out that the collapse of a stochastic collection of dilatonic waves 
(in the Einstein frame) does also lead to an anisotropic pre-big-bang phase . 

In the framework of an anisotropic pre-big-bang phase it is plausible, from a theoretical point of view, to investigate 
the role played by a homogeneous magnetic field which could represent a further source of anisotropy. Thus, the first 
technical point we want to investigate is the possible generalization of the pre-big-bang solutions to the case of a 
homogeneous magnetic field. This generalization is not so straightforward since, in the low energy limit, the dilaton 
field is directly coupled to the kinetic term of the Maxwell field 



R + g ap d a <j>dpct> - ^H^ va H^ a + ^F a0 F al3 



(1.5) 



where F a p = V[ a ^4g] is the Maxwell field strength and is the covariant derivative with respect to the String 



frame metric g^ v . Notice that H^ va is the antisymmetric field strength. In Eq. (1.5) the ellipses stand for a possible 
(non-perturbative) dilaton potential and for the string tension corrections parameterized by a 1 = \ 2 S (notice that 



X s is the string scale). In Eq. (1.5) can be thought as the Maxwell field associated with a U(l) subgroup of 
Eg x E$ . One might think that by going to the Einstein frame the dilaton can be decoupled from the kinetic term of 
the Maxwell fields. This is not the case as we discuss, in greater detail, in Appendix A. 

In addressing the possible occurrence of an anisotropic phase in the life of the Universe there is an immediate 
concern. We want to make sure that the amount of anisotropy encoded in the initial conditions will be eventually 
washed out by the subsequent evolution. It can be shown, in this respect, that an interesting role can be played by 
t he string tension corrections ||. In fact, by adding the first a' correction to the tree-level action reported in Eq. 



(1.5) two interesting things happens. On one hand the curvature invariants get regularized and, on the other hand, an 
anisotropic solution with expanding scale factors gets isotropized. This statement can be made more precise by looking, 
simultaneously at Fig. gand Fig. |j In Fig. ||we plot the shear parameter r(t) = 3[H(t)- F(t)}/[H(t) + 2F (t)} for the 
particular vacuum solution leading to flat axionic spectra in four anisotropic dimensions, namely the case a = —7/9, 
(3 = —4/9. This case does correspond in Fig. [l]to the intersection between the full line and the vacuum ellipse in the 
third quadrant. We can clearly see that for t — > — oo the solution is anisotropic since r{t) is of order one. The crucial 
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question for the purpose of the present paper is by how much the quadratic corrections are able to reduce the degree 
of anisotropy. 

For our future convenience we define the degree of isotropization as the logarithm (in then basis) of the absolute 
value of the shear parameter: 



J(i)= log |r(i) | . 



(1.6) 



The reason for the absolute value in the last equations is that r(t) can change sign depending upon the relative 
magnitude of H and F. In Fig. ^ we plot the degree of isotropization (at the left) and the shear parameter itself (at 
the right). We see that the degree of isotropization is a function of the duration of the string phase. Thus if we have 
a very long string phase we can have I(t) <C —7. On the other hand if the duration of the string phase is not too 
long, then one can naturally have larger values for r(t). Suppose now that the string phase stops at some time t and 
it is replaced by a radiation dominated phase. Then by the continuity of the scale factor s at the transition time the 
anisotropy in the expansion is translated in a temperature ansisotropy according to Eq. (1.1). 




FIG. 2. We plot the evolution of the degree of anisotropy (left) and of the shear parameter itself (right) in the case of 
a mildly anisotropic solution with a = — 7/9 and /3 = —4/9 corresponding to the case of flat axionic spectra |J (i.e. the 
intersection of the thick full line with the vacuum ellipse in the third quadrant of Fig. [I]). As we can see for t — > — oo, 
r M "~ * [~ 7/9 + 4/9]/[— 7/9 — 4/9] = 3/11. For t — > +oo the shear parameter and the degree of anisotropy gets reduced 
as the result of the addition of of higher order string tension corrections (see also Section IV). We notice that the degree of 
isotropization crucially depends upon the duration of the string phase (left). In a short string phase I(t) can be of the order of 
—3, —4. If the string phase is very long the degree of anisotropy can be much smaller than — 7. Notice that, by only looking at 
the right picture we might guess (in a wrong way) that the shear parameter tends (for t — > +oo) to a constant (small) value. 
The crucial point we want to stress is that this is not the case as we can understand by looking at I(t) which is always a 
decreasing function of the cosmic time. 



Therefore, it could happen, depending upon the duration of the string phase, that either I(t)> 10~ 5 or I(t)<10~ 5 
right after the end of the string phase (under the assumption that radiation sets in instantaneously) . Many questions 
can arise at this point. Will the radiation evolution completely erase the primordial anisotropy? If this is not the case, 
should we exclude, on a phenomenological ground, too short string phases since they might potentially conflict with 
the observed value of the large scale anisotropy? Which is the role of a primordial magnetic field in the picture? We 
want to recall, in fact, that even if we would not include a magnetic field in the pre-big-bang phase such a field will 
be anyway dynamically generated by parametric amplification of the electromagnetic (vacuum) fluctuations [§) and 
it will evolve, in the radiation dominated (post-big-bang) phase affecting the evolution of r(i) [fl. Certainly one can 
argue that the magnetic field generated from the electromagnetic vacuum fluctuations is stochastic in nature as all 
the other fields (gravitons, scalar fields...) produced as the result of parametric amplification of quantum mechanical 
initial conditions ||,[l(|. Therefore, we could conjecture that the obtained magnetic field does not have any preferred 
direction. We want to point out that this conclusion has been reached by looking at the magnetic field generation in 
a completely homogeneous and isotropic dilaton-driven phase. If the background is not completely isotropic during 
the dilaton driven evolution, then, we can argue, in analogy with the graviton production in anisotropic backgrounds 
of Bianchi type I, that the produced background has different properties in different spatial directions Jll| . 

Finally, the considerations and the various questions we just stated hold in the case of transition from the string 
phase to the radiation dominated phase. However, it could happen that before reaching the radiation phase the 
background evolves toward a state of decreasing dilaton and then, we should treat again the same issues in the case 
of a background with decreasing dilaton coupling. These are the problems we want to investigate in the second part 
of the present paper. 

Our plan is the following. In Section II we will discuss the magnetic solutions of the low energy beta functions 
whereas in Section III we will discuss their physical interpretation. In Section IV we will investigate the role of the a' 
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corrections and we will try to follow the evolution of the degree of anisotropy all along the string phase. In Section V 
we will study the fate of the anisotropy in a decelerated phase and we will address the issue of possible bounds on the 
duration of the string phase from too large temperature anisotropies. Section VI contains our concluding remarks. 
Some useful technical results concerning the description of anisotropies and of magnetic fields in the String and in 
the Einstein frames are collected in Appendix A. 



II. MAGNETIC SOLUTIONS 



Consider a spatially flat background configuration with vanishing antysimmetric field strangth (H^a = 0) and 
vanishing dilaton potential. The dilaton depends only on time and the metric will be taken fully anisotropic since we 
want to study possible solutions with a homogeneous magnetic field which is expected to break the isotropy of the 
background: 



ds 2 



g^dx^dx" = dt 2 - a 2 {t)dx 2 - b 2 {t)dy 2 - c 2 (t)dz 2 



(2.1) 



This choice of the metric corresponds to a synchronous coordinate system in which the time parameter coincides with 
the usual cosmic time (i.e. goo = 1 and goi = 0). 

By varying the effective action with respect to </>, g^ and the vector potential we get, respectively 



R - g^d^d^ + 2tf#V a Vft 
1 



l 5 u R 

2 M 
— <i> pp-v 



0. 



where V M denote s co variant differentiation with respect to the metric of Eq. (2.1). Inserting Eq. 
we get that Eq. (2.3) can be expressed as 



(2.2) 
(2.3) 
(2.4) 

into Eq. (E3) 



. J? pt/a 



0. 



(2.5) 



Consider now a homogeneous magnetic field directed along the x axis. The generalized Maxwell equations (2.4) and 
the associate d Bi an chi identities can then be trivially solved by th e fie ld strength F yz = —F zy . The resulting system 

as 



of equations (2.2), ( |2.5| ) can than be written, in the metric of Eq. (2.1 

HF + FG + FG + H 2 



2<t> + 2{H + F + G)cj) - 



F' + G z +H + F + G 



B 2 



0, 



= H 2 +F 2 + G 2 



Hcj) = HF + HG + H 2 
Fcp = HF + FG + F 2 - 



H + F + G, 
H, 

B 2 



F - 



2b 2 c 



Gcf) — HF + FG + G + G - 



B 2 



(2.6) 

(2.7) 
(2.8) 

(2.9) 
(2.10) 



where Eq. (2.6) is the dialton equation and Eqs. (2.7)— (2.10) do correspond, respectively to the (00) and (ii) 
components of Eq. (2.5). Notice that B is the magnetic field intensity. The inclusion of a time-dependent B is in 
principle possible but in order to be consistent with the Bianchi indentities we would need to include also an electric 
field (not necessarily homogeneous) . Notice that in the limit B — > we obtain the well known form of the low energy 
beta functions in a fully anisotropic metric. 

In order to give a general solution of the previous system it is convenient to define a generalized "conformal" time 
namely 



irj = dt, 



logV"-5, 



(2-11) 



where <fi is the "shifted" dilaton which is invariant for duality related vacuum solutions. In terms of rj Eqs. (2.2)-( [Oc| ) 
become, respectively, 
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20" + cf>' - 2£' - 20'£ + A 
0" + 0' 2 - 2: 
W' = 0, T' 



1 -Id, 

a e , 



" + 0' 2 - 2£0' + A - £' = 0, 



B 2 



a 2 e~ 2 <t 



where 



S=(W + ^ + g), A = 2(K^ + Wg + ^), 



(2.12) 
(2.13) 
(2.14) 

(2.15) 



with the obvious notation t hat (log a)' = H, (log 5)' = T, (logc)' = Q (the prime denotes differentiation with respect 
to rj). By subtracting Eq. fl2.13| ) from Eq. (fTl|) we obtain 



d2 

E' = -— a 2 e- 2 t 
2 



(2.16) 

Using the remaining equations in Eq. ( 2.16] ) in order to eliminate £' we get a general form of the solution which reads 

a{rj) = a e n ° 7 \ b{ri) = bo^e^, c(r?) = a e e °V, (2.17) 

where satisfies 



2Ho 



A = 0, A = 2(H T + T G + T Go 



(2.18) 



By solving this last equation and by inserting the (trial) solution (2.17) into Eqs. (2.12)-( 2.10D we will get a set of 
consistency relations among the different integration constants which will determine a general form of the solutions 
in terms of the (undetermined) set of initial conditions. 

Let us notice that by repeatedly combining the equations of motion we obtain a very useful relation, namely 
20" = (Ba) 2 exp [—20]. Since the left hand side of this last equation is positive definite, we have also to demand, for 

Ao > 0, which also implies 



physical consistency that, from Eq. (2.18), —0' + 2Ho4>' 



Ho 



H 2 + A < 0' < H - Jn 2 + A 



(2.19) 



Notice that this is a requirement to be satisfied in the presence of a constant magnetic field and it does change in 
the presence of a constant electric field (indeed for a constant electric field directed along the x direction and in the 
absence of any associated magnetic field we would have that, in the same conformal time parameterization 0" < 0). 
The wanted solution is 



a(r?) = a Q e no, \ 6(77) = 6 oe *> e ^ b+WD+A °^ 



e -2A„ V e 4 



c(r?) = Coe ^e (e ° +Wo+Ao)? ' 



e -2A,rn p4 



0(r?) = 0o + (Ho + A )?7 + log + e' 2 ^ 
where, by consistency with all the other equations we have 

Ao = J Hi + ha = 



ao -t 



2V2 



io+0i) ft 



(2.20) 
(2.21) 

(2.22) 



A particular solution of the system is fi xed b y fixing the value of the magnetic field in string units. Before discussing 
the physical properties of the sol ution (2.21) we want to see in which way the vacuum solutions appear in the time 
parameterization defined by Eq. (2.11). The vacuum solutions are quite straightforward using rj and they correspond 
to constant Hubble factors and linear dilaton, namely 



H=H , T = To, G = Go, 0(r/) = O + 2 r/, 



with the condition, 



6 2 = Ho + To + Go ± \ H 2 +T 2 + Gl 



(2.23) 



(2.24) 
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It is straig ht forward to see that these solutions are indeed Kasner-like by working out the cosmic time picture from 



Eq. ( 2.11 ). Indeed, with little effort we can see that 



where 



Ho T Go ,„ „x 

a l = 7, ~Ti a 2 = -p, T - 1 C*3 = — (2.2b) 

Ho + To + Go - 02 " Ho + To + Go - 4>2 Ho + To + Go - 4>2 



Using Eq. (2.24) we can obtain that the thre e (co smic time) exponents satisfy the Kasner-like condition, namely 
a\ + a 2 . + a i — 1 and the ± ambiguity of Eq. ( 2.24 ) simply refers to the two (duality related) branches. 



Before concluding the present Section we want to mention that the equations of motion can also be integrated in 
the case where the coupling of the dilaton to the Maxwell fields is slightly different |Q from the one examined in this 
Section, namely in the case where the action is written as 



S 



R + g ap d a <j>d (i <p + l e -^F af3 F af} 



(2.27) 



with q ^ 0. 

In this case the equations of motion in the generalized conformal time rj can be written as 

20" + 0' 2 - 2S0' + A - 2E' = -i±-lBV e -^ +2 ^ (2.28) 

0" + 0' 2 - 2E0' - £' + A = -|s 2 a 2 e -(«+ 2 ^ (2.29) 

H ' = ! B Y e -( 5 +* p = i±l B 2 a 2 e -( q +2 )( b g , = £±i B a 2 e -fo+ 2 >* (2.30) 
4 4 4 

Notice that in this case the system changes qualitatively but it can still be integrated. In fact, by following the 
procedure outlined in the previous paragraphs we obtain a decoupled equation for (/>' 



g + 1 2( g +l) 2 g + 1 

This equation can be easily integrated providing the solutions to the whole system in the case q =/= 0. 

III. ANALYSIS OF THE MAGNETIC SOLUTIONS 

In order to analyze the previous solution let us assume, for sake of simplicity, that b{rf) = c(i]). Suppose now, that 
our initial state is weakly coupled (i.e. <po sufficiently negative) and (Ho, To) are sufficiently small in string units. 
Even though Ho and To are of the same order they do have different initial conditions. Let us therefore analyze 
what is the role of a magnetic field, sm all in string units. The discussion becomes simpler if we start looking at the 
77 behavior of the solutions. From Eqs. ( 2.2l| ) we can argue that the evolution of 4>'{rj) is quite straightforward since 



for -q —>■ —00 and 0' — > Ho — A whereas for 77 — * +00 we have that 4>'(ri) — * Ho + A . Moreover, as we discussed 
in the previous Section, <fi" > which implies that <f>' is always an increasing function of rj. As a consequence we 
will also have that T(rj) will also be an increasing function of rj whereas H is frozen to its constant value given by 
Ho- According to Eq. ( [2.22 ) the magnetic field intensity in string units determines, together with (f>o the value of 



Ao. Therefore, we conclude that the magnetic field intensity is essentially responsible of the amount of growth of the 
dilaton energy since it turns out that \<f>'(+oo) — <p'(—oo) \ = 2Aq. Thus, by increasing the magnetic field we obser ve an 
increase of the dilaton energy density. Notice, however, that this growth cannot be so large. Indeed, from Eq. ( 2.22| ) 



we notice that the relation between Ao and B is suppressed by the second power of the (initial) coupling constant 
g{4>) = exp(0/2). So, with our initial conditions, Ao is small from the beginning since 0o is very negative. Notice that 
to have small coupling is essentially dictated by the form of our original action which cannot be trusted if g{<j>) gets 
of order one. This qualitative behavior of the solutions is illustrated in Fig. |3|-{| in few specific cases. 
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FIG. 3. We plot the magnetic solutions given in Eqs. (2.21) in the case where TLo = 0.001, B = 01, <f>a = —5. In the left plot 
we illustrate the behavior of TL {dot-dashed line), T (thin line) and (j>' (full thick line). At the right we illustrate the behavior of 
the shear parameter. We notice that, in spite of the appearance, the shear parameter does not decrease, asymptotically but it 
tends to a constant value of order one, in sharp contrast with what happens if we include quadratic corrections to the tree-level 
action (see next Section). We see that TL, J- and 4>' are almost constant if the magnetic field is small in String units. 

As a a reference value for TLo and for <po we take 0.001 and —5 since we want to deal with small curvatures in String 
units and small couplings. In Fig. || (left) we report the solutions in the case of a quite large magnetic field B ~ 1. 
The dot-dashed line denotes the behavior of TL (constant), the full (thin) line denotes the evolution of F(r]) whereas 
the behavior of the full (thick) line denotes the evolution of (/>' . 
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FIG. 4. We plot the same solution given in Fig. |3| but with a different value of the magnetic field which is now B = 0.001. 
As in the previous plot at the left we report the evolution of TL (dot-dashed line), J- (thin /ine)and (f>' ( full thick line). Notice 
that in this case the growth of <f)' is so minute that the dot-dashed line and the full thick lines are almost indistinguishable. 
This is due to the fact that 4>' continuously interpolates between TLo — Ao and TLo + Ao and Ao ~ O.OOle -5 . At the right we 
report the evolution of the shear parameter for the same solution. 

It is interesting to compare this plot with the ones illustrated in Fig. |^ where (left plot) we report the solution 
(2.21) with the very same choice of initial conditions of Fig. (|J) but with different magnetic field, namely B = 0.001 
in String units. We see that the growth of </>' is reduced (by two order of magnitude). 

What about the evolution of the shear parameter? Let us look at the left plots reported in Fig. |^, ^ where the 
evolution of 



r(r]) 



(3.1) 



is given for the same set of parameters we just discussed. In Fig. || (where B ~ 0.1) we would say that the shear 
parameter decreases. This is certainly true but if we look at the numbers we can see that r(rf) starts of order one 
and essentially remains of order one. In Fig. |] (i.e., B = 0.01) r(rf) seem to decrease, but again, a more correct 
statement would be that r(rf) remains of order one (i.e. r(+oo) ~ —0.88), up to a transition period. Thus, in spite 
of its smallness, the magnetic field has always the property of conserving the anisotropy. The negative sign in r(rj) is 
only due to the fact that T > TL for any r\. Thus, as far as the tree-level evolution is concerned we can say that no 
isotropization is observed as a consequence of the inclusion of the magnetic field. This situation should be contrasted 
with what happens when the string tension corrections are included. 

The second effect we see is that the smaller is the value of the magnetic field, the milder is the transition to the 
vacuum solutions. We said, in the previous section that the vacuum solutions are essentially given by TL — constant, 
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T = constant, <f>' — constant. For n — > +00 the magnetic solutions should then tend to the vacuum solutions. Thus, 
the influence of the magnetic field is mainly in the transition period from the initial state to the (asymptotic vacuum) 
solution. 

In order to show that the large n behavior of our solutions is Kasner-like, let us look explicitly to the cosmic time 
description of the system. Let us l ook a t the evol ution of our system in cosmic time. The relation between n and t 
can be obtained by inserting Eqs. ( 2.21 ) into Eq. ( 2.11 ) whose integrated version becomes 



t 



ao&§e*> 



[2{T + Wo) + A ] [2(Jb+Wo)-A ]' 



(3.2) 



Notice that the only two indep ende nt quantities determining the relation between n and t are H.q and Ao (which is 
related to B according to Eq. ( 2.22 )) since To can be expressed in terms of Ao and Jio- 

If we now take the large 77 behavior of Eq. ( |3.2| ) and we substitute back into Eqs. (2.21) we do find that the two 
scale factors can be expressed, in cosmic time, as 



o(i)~t a , b{t)~lf 



where 



a 



Ho 



2(H + To) + A ' 







F0+H0 + A 
2(H Q + To) + A ' 



(3.3) 



(3.4) 



We can notice that a 2 + 2/3 2 = 1. This can be ea sily seen by recalling the definition of Aq = Ti 2 + 2(2HqTo + Tq) 
which can be deduced from Eqs. (2.22) and ( [2.18 ). 

So, close to the singularity, for n — > 00 the solutions are Kasner-like. Notice that the presence of the vector field 
(with the specific four-dimensional coupling required by the low energy effective action) does not lead to any type of 
oscillatory stage. This is evolution is in sharp contrast with one obtains in the case of a vector field (included following 
the Kaluza-Klein ansatz) in the Einstein frame |13|. The reason for this difference stems from from the difference in 
the couplings. In Fig. || the evolution of (j>, H(t) and F(t) are illustrated in the case of B ~ 0.1. Notice that t — > 
does correspond to the large n behavior. The singularity is located, in this example for t ~ 2. 





FIG. 5. We illustrate the behavior of H(t), <j> (left plot) and of F(t) (right plot) in the case of a magnetic field of the order 
of B ~ 0.1 in String units. 

In Fig. H we report the same solution illustrated in Fig. |j] but with a smaller magnetic field, B = 0.01. As one can 
expect the singularity and the vacuum limit is reached before than in the previous case. 
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FIG. 6. We illustrate the behavior of H(t), </> (left plot) and of F(t) (right plot) in the case of a magnetic field of the order 
of B ~ 0.1 in String units. 



IV. MAGNETIC SOLUTIONS AND THE HIGH CURVATURE REGIME 



In the previous two Sections we examined fully anisotropic (magnetic) cosmologies using the tree level action. In 
this Section we are going to extend and complement our results with the addition to the tree-level action of the first 
string tension correction. In the presence of the first a' correction the action becomes 



R + Q ap d c 



4 



R, 



GB 



(4.1) 



where Rq B is the Gauss-Bonnet invariant expressed in terms of the Riemann, Ricci and scalar curvature invariants 



Rgb — Rfj.vapR^'"* 13 — ^R^yR^ + R 2 , 



(4.2) 



and a; is a numerical constant of order 1 which can be precisely computed depending upon the specific theory we 
deal with (for instance lo = —1/8 for heterotic strings). In principle the magnetic field should also appear in the 
corrections to the tree-level action. However, it turns out that, the terms which would be significant for the type of 
configurations examined in this paper, the gauge fields appear beyond the first a' correction. For example terms like 
pa& piiv n a p^ v w jii appear to higher order in a 1 . Moreover, possible terms like F al3 R a are vanishing in the case of 
constant and homogeneous magnetic field directed, say, along the x direction. Finally, other possibly relevant terms 
(involving contractions of the vector potentials with the Ricci or Riemann tensors) are forbidden by gauge invariance. 

Of course, the spirit of our analysis of the high curvature regime is only semi-quantitative. After all, the first 
curvature correction can be thought as illustrative, since, ultimately all the a' corrections will turn out to be important. 
The only way of turning-off the possible effect of curvature corrections is to have explicit solutions whose curvature 
invariants are already regularized at tree level like in the examples presented in [ [l4| where it was shown that there 
exist weakly inhomogeneous solutions of the tree level action which are regular and geodesically complete without the 
addition of any curvature or loop corrections to the tree level action. If this is not the case, the curvature corrections 
have to be certainly included. The inclusion of higher order curvature corrections can be studied either in the Einstein 
frame jl5| or in the String frame (l6) . Another interesting approach has been outlined in jlTj . 

For a reasonably swift derivation of the equations of motion it is useful to write the action in terms of the relevant 
degrees of freedom. In order to get correctly the constraint equation (i.e. the (00) component of the beta functions) 
it is appropriate to keep the lapse function N(t) it its general form 



. 9AlI , = diag[iV(t) 2 ,-e 2Q « 



,2/3(t) 



2/3(t) 



(4.3) 



where we parameterized the two scale factors with an exponential notation. Only after the equatio ns o f motion have 
been derived we will set N(t) — 1 corresponding to the synchronous time gauge. In the metric (4.3) the previous 
action (4.1) becomes, after integration by parts, 



S = 



1 

2A 2 



dte a+2 ^ 



1 



-0 2 - 2/3 2 - 4d/3 + 2d0 + 4/30 - -B 2 e' 4m 



4/V 3 



d /3 2 - 4 



(4.4) 



By varying Eq. (4.4) with respect to the lapse function N(t) and imposing, afterwards, the cosmic time gauge we get 
the constraint 



^ + 2/3 2 + 4d/3 - 2d0 - 4/30 - 



1 



-B l e 



-4/3 _ 



6do 



0. 



(4.5) 



where we took string units A s = 1 and u> = 1. By varying the action with respect to a, (3 and we get, for N(t) = 1 
the diagonal components of the beta functions and the dilaton equation, 



4/3- 20+(4/3- 20)(d + 2/3 



(a + 20- 0)0/3 2 + 0/3 2 + 2$(3<j) 



+ L(t) = 0, 



4(/3 + a - 0) + 2B 2 e' if) + 4(d + /3 - 0)(d + 2/3 - 



/3d0 + /3d0 + /3d0 + /3d0(d + 2/3 - 



d - 2/3) + 2(0 - d - 2/3) (d + 2/3 - 



2f3(af3 + 2a(3) + (2d/3 2 - 3 )(d + 2/3 - 0) - 30 2 



(4.6) 

+ 2L(t) = 0, (4.7) 
-L(t) = 0, (4.8) 
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where we defined 

L(t) = -0 2 - 2/3 2 - 4d/3 + 2a<j) + 4/30 - ^B 2 e~ Af3 + i(8^d/3 2 - 4 ). (4.9) 

We can numerically integrate this system. The technique is very similar to the one used in the case of fully anisotropic 
string cosmologies with zero magnetic field [p|Jlr|. Our results are illustrated in Fig. ||and[9| There are two physically 
different cases. 
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FIG. 7. We report the result of the numerical integration of Eqs. ( |4.8| ) in the case where the initial conditions are mildly 
anisotropic. We choose B = 0.01. 



We integrate this system by imposing, as initial conditions small curvature and small dilaton coupling. We also 
impose mildly anisotropic initial conditions (corresponding to a > 0, b > and a > 0, b > 0). One of the results 
stressed in the previous Section is that the tree-level solutions with magnetic field are singular. Moreover, close to the 
singularity no oscillatory behavior is present. Therefore, as time goes by and as the solutions approach the singularity 
the vacuum solutions are recovered. One can then argue that if the a' corrections become relevant when the vacuum 
regime has been already recovered, nothing should change with respect to the case where B = 0. This in some sense 
is what happens, but not exactly. 



2.25 




FIG. 8. We illustrate the numerical solutions of the equations of motions with string tension corrections and we plot the 
shear parameter. We select expanding initial conditions (in the String frame). In this particular example (left plot) we choose 
the magnetic field to be B = (dot dashed line), and B — 0.01 (full line). The difference between the two cases is quite minute 
and in order to show it we do not plot the full evolution of r(t) (which goes to zero for t > 0). At the right we choose B = 0.01 
and we change the initial condition within the vacuum solutions. The dot dashed line is for b(t) ~ i _4//9 and the full line for 
b(t) ~ t~ ' 2 . Also in this second case the difference is quite minute. We decided to plot only a limited range of time steps in 
order to stress the numerical difference between different values of the initial shear and of the primordial magnetic field. The 
full picture looks like the ones reported in Fig. bl 



If the magnetic field is large in string units (i.e. B > 1) this system evolves towards a singularity, as expected. In 
this case, it is certainly true that the tree level solutions evolve towards their vacuum limit. However, in this limit 
the solutions get more and more anisotropic. So the hypothesis B > 1 simply contradicts the assumption of mild 
anisotropy. From our point of view it is more interesting the case with B < 1. In this case mildly an isot ropic initial 
conditions will be attracted towards an isotropic fixed point. One can understand this looking at Eqs. ( |4.8| ). For t — > 
b(t) is monotonically increasing (if we select, a s we do, expanding initial conditions). Now the magnetic field is always 
suppressed by b~ 4 (t) = e -4 ^ 4 - 1 terms in Eqs. (4.S). Therefore, an initially small magnetic field will become more and 
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more sub-leading. The solutions will then be attracted towards the quasi-isotropic fixed point a ~ f3 ~ 0.606 and 
- 1.414. 




FIG. 9. We report the evolution of the shear parameter in the String phase. With the dot-dashed line we have the case 
B = 0, with the full line the case B — 0.01. If a magnetic field is present the shear parameter is larger than in the case of zero 
magnetic field. It is of crucial importance to notice, for our purposes, that the shear parameter does not tend towards a small 
(finite value) but it decreases towards 0. In order to stress this aspect we also plot the logarithm (in ten basis) of the modulus 
of the shear parameter. We can see, as previously discussed, that the duration of the string phase is proportional to the drop 
in r(t). 

In this case the curvature invariants are finite everywhere. We say that the fixed point is quasi-isotropic since the 
amount of anisotropy depends upon the duration of the high curvature scale. For example, as illustrated in Fig. |[ it 
can well happen that for a very short stringy phase the amount of anisotropy measured by the shear parameter r(t) 
will be of the order of 10 -3 . 



V. THE FATE OF THE ANISOTROPIES IN THE POST-BIG-BANG EVOLUTION 



In this Section we would like to investigate the fate of the anisotropics possibly present in the dilaton-driven and 
in the String phase. In particular we would like to understand in which limit these anisotropies can be either reduced 
to an acceptable value or completely washed out. For numerical purposes we find useful to exploit, in the present 
Section, the Einstein fra me p ic ture . The evolu tion equations in the Einstein frame are obtained in Appendix A. By 
linearly combining Eqs. flA.gj ), ( |A.9| ) and ( A. 10 ) of Appendix A we obtain a more readable form of the system, namely 



H- 

04 



H(H + 2F) 



P 
6 



w 



. F {H + 2F) = P - + ^e 
(H + 2F)4> = we-*, 

+ 2Ff - (H 2 + 2F 2 ) = ^- 
-(H + 2F)p = 0, to 



2 

4Fw 



0, 



w 



26 4 ' 



(5.1) 
(5.2) 
(5.3) 
(5.4) 

(5.5) 



where we assumed the presence of a radiation fluid p = p/3 accounting for the field modes excited, via gravitational 
instability, during the dilaton driven phase Jl§| ] and re-entering in the post-big-bang phase. It is in fact well known 
that the ultraviolet modes of a field parametrically amplified behave as a radiation fluid |ll| . 

From Eqs. (|5.l|)-(5.5) we can directly obtain the evolution equations for the quantities we are interested in, namely 



"(*) = 



3(H - F) 
H + 2F ' 



n{t) 



H + 2F 



g(<) 



A*) 
p(t)' 



(5.6) 



where n(t) is the mean expansion parameter r(t) is the shear anisotropy parameter and q(t) the fraction of magnetic 
energy in units of radiation energy. Notice that q(t) does not correspond to the critical fraction of magnetic energy 
dens ity since, in principle, we have to take into acc ount the energy density of th e di laton field. By subtracting Eq. 
(5.2) from Eq. (O) and by using the constraint ( |5.4| ) together with the definitions ( |5.6| ) we get the shear evolution By 
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summing Eqs. (5.1) an d (|5.2| ) with similar manipulations we get the evolution of n. The resulting system, equivalent 
to the one of Eqs. (5.1)-( |5T5f ), but directly expressed in terms of n, r and q reads 



nr + rn + 3n r = —qpe 



3n 



9?i 2 = |(l + ge- 



-nrq 



6n 2 (l 



-) = ^-+p(l + «fe~*), p + 4np = 0, 4> + 3n<}> = pqe-*. 



(5.7) 
(5.8) 

(5.9) 



We want now to study the combined action of the dilaton and of a radiation fluid in a post-big-bang phase. Let 
us start from the case of medium size anisotropies. Suppose, in other words that after a string phase of intermediate 
duration the shear parameter is of the order of 10~ 3 . As we saw from the previous Section this value is not excluded. 
Before analyzin g th e gene ral case let us recall the results of the case where the dilaton is absent. In this case the 
system of Eqs. (5.7)— (5.9) can be simplified as follows 



nr + rn + in r 



-qp, 



p 4 
3h + 9n 2 = -(1 + q), q — -nrq = 0, 

r 2 

6n 2 (l - — ) = p(l + <7e~*), p + Anp = 0. 



(5.10) 
(5.11) 

(5.12) 



We are interested in the case of magnetic fields which are under-crtical , namely q(to) < 1 where to is the initial 
integration time. Suppose, moreover, that r(to) ~ 10~ 3 . Needless to say that this is just an illustrative value for r 
and it is not meant to be of any particular theoretical relevance. As we stated clearly in the previous Section (see 
for instance Fig. || ) if the string phase is long the initial shear parameter can be much smaller. The results of the 
numerical integration are summarized in Fig. [l^ where we report the logarithm (in ten basis) of the modulus of the 
shear parameter. 
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FIG. 10. We report the results of the numerical integration of Eqs. ( 5.10 )— ( 5.12 ). We take r(to) — 10~ and we integrate the 
system forward in time for different values of q{to)- In the left plot we illustrate the evolution of the logarithm (in ten basis) of 
the modulus of the shear parameter. With the dot-dashed (thick) line we represent the case q — where the shear parameter 
is known to decrease sharply as t~ 1//2 . If, initially, the critical balance between the magnetic energy densit y and the radiation 
energy density increases, then r gets attracted towards and asymptotic value as described by Eqs. ( |5.13| ). At the right (full 
thick line ) we report the evolution of q for the case q(to) = 0. 01. With th e thin line we report the qualitative estimate based 
on Eqs. (5.13) and obtained by solving, approximately, Eqs. (p.l0[)-(p.l2[) for r < 1 and q < 1. 



At the bottom of the left plot, with the dot-dashed line, is reported the case o f zero magn etic field (i.e. q(to) = 0). 
This is the simp lest case since we can approximately solve the system of Eqs. ( 5.1C )— (5.12) in the limit of small r. 
From Eq. ( 5.10f ) we can find that a consistent solution is 



r 
~2t 



3q 

t ' 



rq. 



(5.13) 



So, if q = we can clearly see that r(t) ~ t -1 / 2 . This is nothing but the well known result that in a radiation 
dominated phase the shear parameter decreases as l/y/i and it is precisely what we find, numerically, in the dot- 
dashed line of Fig. [To] (left plot). If we switch-on the magnetic field we also know, from Eqs. (5.13) that the anisotropy 
will not decrease forever but it will reach an asymptotic value which crucially depends upon the balance between the 
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magnetic energy density and the radiation energy density. In fact from eqs. ( |5.13 ) we see that r — for \r(t) \ — > 6|q (t)| 
and this is precisely what we observe in the full lines of Fig. |h] where the integration of the Eqs. ( 5.10| )— ( 5.12 ) is 
reported for different (initial) values of q. We see that the asymptotic value of attraction for r is roughly six times 
the in itial value of q. The evolution of q(t) when f(t) — ► can be simply obtained by integrating once the second of 
Eqs. ( |5.13D for the case r ~ — 6q. The result is that, taking for instance q(t Q ) = 0.01, q(t) ~ l/{41og[(i/to)] + 100}. 
This simple curve is reported in Fig. |l^ at the right (full thin line). As we can see the full thick line (result of the 
numerical integrations is practically indistinguishable). 

The conclusion we draw for our specific case is very simple. If the string phase is not too long and if sizable 
anisotropies (i.e. r ~ 10~ 3 ) are still present at the end of the stringy phase, then , in the approximation of a sudden 
"freeze-out" of the dilaton coupling, the left over of the primordial anisotropy will be washed out provided the 
magnetic field is completely absent. If a tiny magnetic seed directed along the same direction of the anisotropy is 
present, then in spite of the primordial anisotropy, the shear parameter will be attracted to a constant value essentially 
fixed by the size of the seed. 

Let us now investigate the opposite case. Let us imagine of switching-off completely the radiation background. 
Then the only non trivial variables determining the shear para m eter and the mean expansion will be the dilaton 
coupling and the magnetic energy density. In this case Eqs. (5.7)— (5.9) do simplify as follows 

hr + rn + 3n 2 r = — we~^, (5-14) 

3n + 9n 2 = ^e~ , w + 4(n - ™)w = 0, (5.15) 



w + tt)w = 0, 

we~^, <j> + 3?i0 = pqe~ 



9 ' 2 

As in the previous case let us assume, in order to fix our ideas and in order to have a swift comparison 
10 -3 . The results of our analysis are reported in Fig. jO , 



(5.16) 
that 




FIG. 11. We report the resluts of the numerical integration of Eqs. ( 5.14 )— ( 5.16 ) for different initial values of the magnetic 
energy density in Planck units. As we can see (dot-dashed line) the shear parameter remains fixed to its initial value if the 
magnetic field is completely absent. As soon as the magnetic field increases the shear parameter gets attracted towards a fixed 
point whose typical anisotropy is of order one. As in the previous cases we took r(io) ~ 1CP 3 and we also took 4>(to) ~ —0.1. 



Suppose first of all that the magnetic field is switched-off, then, as we see from the dot-dashed line the anisotropy 
is completely conserved. The introduction of a tiny magnetic seed will eventually make the situation even worse in 
the sense that the anisotropy will eventually grow to a constant value which depends on the magnetic field intensity. 
In Fig. [ll] we see (full lines) that already with a magnetic energy density of the order of 10 -3 (in Planck units) the 
increase in the anisotropy cannot be neglected. Therefore, if the transition from the string phase to the radiation 
dominated phase does occur through an intermediate dilaton dominated phase of decreasing coupling we have to 
accept that the anisotropy will be conserved provided the magnetic field is absent. If a magnetic field is present the 
anisotropy gets gets frozen to a constant value which can be non negligible depending upon the size of the magnetic 
seed. In very rough terms our analysis seems to disfavor mildly anisotropic models with a short string phase and 
with the simultaneous occurrence of a long dilaton dominated phase prior to the onset of the radiation epoch. So 
the combination of short string phase, long dilaton driven phase and radiation phase shorter than usual might give 
too large shear, say, prior to nucleosynthesis. Our conclusion can be even more problematic if it is present a sizable 
magnetic field. 

W it h th e knowledge coming from the two previous examples we can easily analyze the general case given in Eqs. 
( |5.7| )-(5.S). The dynamics of the system will depend upon the different balance in the initial conditions. Let us 
suppose, firstly that the dilaton kinetic energy is slightly smaller (by two orders of magnitude) than the radiation 
energy density. The results of the numerical integration of the shear parameter are reported in Fig. 12. 
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FIG. 12. In this two plots we report the results of the integration of Eqs. ( 5.7 )—( 5.£ ) in the general case where magnetic field, 
dilaton and radiation fluid are present simultaneously. In these particular plot we assume that, initially, the dilaton energy 
density is smaller (by two orders of magnitude) than the radiation energy density. The left plot corresponds to the case of 
4>(to) = —0.2 and the right plot corresponds to the case of cj>(to) = —2. Different initial q are reported in both cases. As in the 
previous cases we assumed r(to) ~ 10 -3 . 



We can clearly see that by increasing the magnetic field the anisotropy will increase. The dilaton evolution clearly 
affects this process as we can argue by comparing Fig. |lj with Fig. [l(] where the dialton was absent. It is also 
interesting to notice that the initial value of the coupling can be of some relevance. In Fig. |l2] we reported the result 
of the integration for two values of the i nitia l coup ling, namely (f>o = —0.1 (at the left) and 4>q — —2 (at the right). 
The effect can be explained from Eqs. ( 5.7 )— ( 5.9 ) where we can see that by tuning the initial coupling to smaller 
values we amplify the effect of the magnetic energy density which appears always as qe~^. 

If the radiation energy density is sub-leading with respect to the dilaton kinetic energy we can expect, on the basis 
of the intuition developed in the previous cases that the shear parameter will stay basically constant if the magnetic 
field is turned off and it will increase to a constant value. This is more or less what happens. In Fig. [l3] we report 
the results of a numerical integration of the shear parameter in the case where the radiation energy is hundred times 
smaller than the dilaton kinetic energy. 




FIG. 13. We report the numerical integration of the system given in Eqs. (^^)-(^9|) for another set of parameters. In this 
case we assume that the radiation energy density is hundred times smaller than the dilaton energy density and we watch the 
relaxation of the anisotropy for different values of the initial magnetic fields. As in the previous case the left picture refers to 
the case where (j)(to) ~ —0.2 whereas the right plot refers to the case where c/>(to) ~ —2. These plots have to be compared with 
the ones reported in Fig. |l(] and Fig. [ll]. We see that, unlike in the case of Fig. [Hj, the integration with q(to) = does not 
lead to an exactly constant anisotropy. For larger values of q(io) the shear parameter saturates, as expected. 



The left picture refers to the case of <f>(to) ~ —0.2 and the right figure refers to the case of <j>(to) ~ —2. In both 
cases we can notice that in the limit of zero magnetic field the anisotropy is not exactly constant as in the case of 
p = reported in Fig. [ll] but it has some mild slope and it decreases. At the same time, in contrast with Fig. |l^ the 
asymptotic value of the anisotropy is not exactly the one we could guess on the basis of the argument illustrated in 



Fig. 10 



We want to stress, finally an important point. Suppose that the dilaton kinetic energy and the radiation energy 
are both present but the magnetic field is zero (q(to) = 0). Then we said that the anisotropy has a mild slope. The 
question is how mild. This point is addressed in Fig. ^ where the results of a numerical integration are reported for 
fixed dilaton kinetic energy and for radiation energy ten, hundred, thousand ad ten thousand smaller than the dilaton 
kinetic energy. 
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FIG. 14. In this plots we assume that the magnetic field is initially zero. We also assume that the radiation energy density 
is smaller than the dilaton energy density which is assumed to be of order one. The plot at the left is for </>(io) = —2 whereas 
the plot at the right is for the case 4>(to) = —10. Nothing changes in the two plots since, for both cases q = 0. We notice that 
the behavior of the anisotropy is qualitatively similar both to the behaviors reported in Fig. [To] and Fig. ^ but with crucial 
differences. For instance, take the case with p ~ 0.1 at the left. In this case the radiation and dilaton energy densities are 
almost comparable. So there is some functional similarity with the plot of Fig. |^. Look anyway at the time scale and take 
into account that the initial value of the anisotropy is the same in both cases. The same reduction of the shear parameter 
obtained in 100 time steps in the case of pure radiation is now obtained in 10 6 time steps. So the dilaton slows down the decay 
of the anisotropy in the absence of magnetic field. If the radiation energy density is gradually removed from the system the 
anisotropy is ( almost) constant. 

We see that, in the best case, the anisotropy decrease of two orders of magnitude in 10 6 time steps. This behavior 
should be compared with the sudden decrease experienced by the anisotropy in the case of Fig. |l^ (dot-dashed 
line) where the dilaton was absent. Therefore, the presence of the dilaton in a radiation background slows down 
considerably, in the absence of magnetic field, the sudden fall-off of the shear parameter. 

The conclusion we can draw from the analysis of the general case is that, also in the presence of the dilaton field, 
the magnetic field has the effect of increasing the anisotropy. Moreover, if the dilaton is larger than (or comparable 
with) the radiation fluid at the moment of the end of the string phase the anisotropy does decrease not so easily. 

VI. CONCLUDING REMARKS 

In this paper generalized the solutions of the low energy beta functions to the case where a magnetic field is present. 
We found that these solutions can be expressed analytically in the String frame. We investigated their limit in the 
vicinity of the singularity and we noticed that they approach, monotonically, the well known Kasner-like vacuum 
solutions. 

We addressed a very simple question which can be phrased in the following way. Suppose that the dilaton-driven 
evolution of the string cosmological models is fully anisotropic thanks, for instance, to a primordial magnetic seed. 
The obvious issue is to understand how (and possibly when) these fully anisotropic solutions reach the complete 
isotropic stage. By complete isotropic stage we meant, more quantitatively, a value of the shear parameter smaller, 
say, than 10~ 7 . We reached a number of conclusions. 

For physical reasons we want to deal, from the very beginning, with anisotropic solutions whose scale factors are 
both expanding and accelerated. In this situation it is not forbidden to have large anisotropies and shear parameters 
of order one prior to the onset of the string phase. 

Now, if the string phase is sufficiently long and if it is immediately followed by a radiation dominated phase, then, any 
pre-existing shear coming from the initial condition is very efficiently washed out to arbitrary small values depending 
upon the duration of the string phase. Of course, this conclusion holds provided no magnetic field is present in the 
ordinary decelerated phase. If a magnetic field is present, then the amount of shear appearing in our present Universe 
will not be determined by the primordial shear but by the shear induced by the magnetic seed itself. The role of the 
magnetic field might also be of some help but only in the case where the primordial shear is really large after the 
string phase. In this case a very small magnetic field can effectively attract the large anisotropy towards a smaller 
value. 

We can also have the opposite situation. Namely we can have the case where the string phase is very short and 
it is followed by a phase dominated by the dilaton until sufficiently small scales lasting from the end of the sting 
phase until the onset of radiation which should anyway occur before nucleosynthesis. In this second scenario the shear 
will not be erased for two reasons. First of all because the duration of the string phase is proportional to the shear 
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reduction. Secondly, because in a dilaton dominated phase (in the absence of magnetic field) the primordial shear is 
conserved. 
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APPENDIX A: FROM THE STRING TO THE EINSTEIN FRAME 



In this Appendix we show explicitly how to get the Einstein frame action for our system. The transformation from 
the String to the Einstein frame does not remove the coupling of the dilaton to the kinetic term of the gauge fields. 
The String frame action is simply 



S 



d x\J—ge 



(A.l) 



In four dimensions the transformation from the String frame metric (g^v) to the Einstein frame metric (G^ v ) simply 
reads 



9^v 



-G. 



Therefore, 



R = e" 



-ge 



K-3G af3 V a V^ 
b F rF G ap <fP - 

r af3 r pcrg g — 



^Ge-^F^F^G^G^, 



(A.2) 

(A.3) 
(A.4) 



the derivatives in Eq. ( |A.3| ) are covariant with respect to the metric G^ v . Thus, the Einstein frame action can be 
written as 



S F = 



d i x\^G 



-U + -G af3 d a <j)d <t) - -e-^FaaF^ 



(A.5) 



Again from this equation we can derive the Equati ons of motion. We will just report the essential points. The 
Equations of motion derived from the action of Eq. (A.5) read 



^ 2 2 



TZ(d) + e-*T»(m)+T»(r) 



*F n «F a P 



-(ppafl 



(A.6) 



with 



l -5^d a 4>dp4>, r;(m) = -f^f™ + -^f^f^ , t;(/) = dia g ( P , - P , - P , - P ), (A.7) 



where on top of the energy momentum tensors of the dilaton (i.e. T^(d)) and of the Maxwell fields (i.e. T^(m)) we 
also added the energy momentum tensor of the fluid s ourc es (i.e. T"(f)) because of the considerations reported in 
Section V. In the fully anisotropic metric given in Eq. (2.1) the previous equations of motion become 



3F 2 + 2F 



4+ 










4 




B 2 








3' 


+ F 2 


f HF = 


4 



p_ B*_ 
2 ~ 46 4 ' 



(H + 2F)0 



B 2 



2 6 4 ' 
p+ (H + 2F)(p + p) = 0. 



(A.8) 

(A.9) 

(A.10) 

(All) 
(A.12) 



Concerning these equations two technical comments are in order. First of all the over-dot denotes the derivation 
with respect to the cosmic time of the Einstein frame which is related to the cosmic time of the String frame (used, 
for instance in Section II) as 



dt E = e ^dts, a E (t E ) = e *a s (t s ), b E (t E ) = e H s (t s ). 



(A.13) 
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In the same way H E = (log a)' and F E — (log b)' are the Hubble factors defined in the Einstein frame. We will not 
denote explicitly this distinction but we would like to remind that the various quantities defined in Sections II, II, 
and IV are defined in the String frame whereas the discussion reported in Section V refers to the Einstein frame. 

Most of the times one can check that all the "physical" quantities are invariant with respect to a change of frame. 
Typical examples of this property are the spectra of the fields excited by the dilaton growth in the case of the Kasncr- 
like "vacuum" solutions derived in Section II. Suppose in fact to compute the axionic (or gravitonic or dilatonic) spectra 
in the String frame. Then Suppose to redo the same exercise in the Einstein frame with the transformed solutions. 
The two spectra will be equal. It is curious to notice that the shear parameter is only approximately invariant under 
conf orma l rescaling. In order to show this aspect let us consider the dilaton-driven (vacuum) solutions reported in 
Eq. (1.4) and written in the String frame 

+ 1/3 



a s (t s ) = 



b(t s ) 



</>(t,) = (a + 2/3-1) log 



(A.14) 



with a 2 + 2f3 2 = 1. The shear parameter for this solution can be very simply computed in the string frame 

3[H s (ts) - Fs(t s )] = 3(a~(3) 
n s) [H s (t s ) + 2F s (t s )} ~ a + 2p ■ 

Let us therefore do the same exercise in the Einstein frame. The scale factors and the cosmic time, transformed from 
the String to the Einstein frame are 



(A.15) 



a E {t E ) 
Therefore we get that 



t_E 



a-20 + 1 
3-a-2p) 



b E {t E ) 



t E 


3-a-2£ 


t. 




t E 


3-a-2p 


*1. 


) 


h. 




h. 





r{t E ) 



3[H E (t E ) - F E (t E )] _. 6(a - 0) 
[H E (t E ) + 2F E (tE)} 



(A.16) 



(A.17) 



3 - (a + 2/3) 

Now, by comparing Eqs. ( A.15 ) and ( A.17 ) it is clear that r s and r E go to zero in the same way and they are both 
proportional to a — j3. In this sense we can say that r is approximately invariant under conformal rescaling. For 
example by taking the anisotropic solution leading to flat axionic spectrum namely a — —7/9, [3 — —4/9 we 
get that |r s (i s )| ~ 3/5 whereas |rE(tg)| ~ 3/7. In all the cases discussed in this paper we explicitly checked that 
the shear parameters in one frame and in the other are of the same order of magnitude. So for practical purposes 
the shear parameters computed in the cither in the String or in the Einstein frame have the same quantitative 
information provided the calculation is performed with respect to the correct cosmic time, which changes from frame 
to frame. More generally one would like to find a frame-independent measure of the degree of anisotropy. One could 
argue that these quantity is provided by the Weyl tensor and, in particular by C J* ' . Indeed it is easy to show that 
C = Cftvof where the calligraphic style denotes, as usual in this Appendix, the quantities computed in the Einstein 
frame. Notice that the index position is crucial j2Q] since C^yap — e^C^ap. The only problem with G„ VC P is that it 
is dimension-full. If we want to construct some dimension-less combination we run into the same ambiguity we just 
described since C^afsC^ ^ 13 / R flua fjR muual3 transforms non trivially under conformal rescaling. So our conclusion is 
that the shear parameter is perhaps still the best quantity in order to characterize the degree of anisotropy. In our 
problem, moreover, we can show that the way the shear parameter goes to zero is a truly frame independent statement 
since, ultimately if r s (t s ) <x \a — /3| we also have that r E (t E ) °c \ot — /3|. 

The second point we would like to stress concerns the covariant conservation of the (total) energy-momentum tensor 
in the Einstein frame. The Bianchi identities impose 



V„ 



T;(<0 + e-*2£(m) + !£(/) 



0. 



(A.18) 



Now, it is easy to see that the component of this equation, in the case of a magnetic field directed along the x axis, 
can be written as 

(9 T °(m) + 4FT °(m))e-^ + U + <p 2 (H + 2F) - T °(m) e -^ + p + (H + 2F)(p + p) = 0. (A.19) 

cancels. Thus, in order to satisfy the 



Using now Eq. ( |A.ll ) we see immediately that the mixed term Tq (m)e 
covariant conservation of the energy-momentum tensor we have to impose 

d T$(rn)+4FT$(m) - 0, p + (H + 2F){p +p) = 0, (A.20) 

which implies that Tq = B 2 /(2b A ) (where B is a constant) as correctly reported in Eqs. ( A.S)-( A.12| ) from the very 
beginning as a consequence of the solutions of the equation for the field strength. 
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